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Abstract. We review some partial results for two strictly related problems. The first problem consists 
in finding the optimal joint unitary transformation on system and ancilla which is the most efficient in 
programming any desired channel on the system by changing the state of the ancilla. In this respect we 
present a solution for dim(J^) = 2 for both system and ancilla. The second problem consists in finding 
the optimal universal programmable detector, namely a device that can be tuned to perform any desired 
measurement on a given quantum system, by changing the state of an ancilla. With a finite dimension d 
for the ancilla only approximate universal programmability is possible, with d — increasing function 

versus e _1 . We show that one can achieve d(e^ 1 ) polynomial, and even linear in specific cases. 
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1 Introduction 

1 A fundamental problem in quantum computing 
and, more generally, in quantum information process- 
ing pQ is to experimentally achieve any theoretically de- 
signed quantum channel or detector using a fixed device, 
through a suitable program encoded in the state of an 
ancillary system. While a large branch of theoretical re- 
search in quantum information addressed the design of 
algorithms and of circuits to solve precise problems, a 
parallel research line is that of designing devices that can 
be programmed to achieve different tasks, just like classi- 
cal computers do. Moreover, designing a programmable 
quantum gate or detector is a problem of relevance for 
example in proving the equivalence of cryptographic pro- 
tocols, e. g. proving the equivalence between a multi- 
round and a single-round quantum bit commitment 0, 
or when trying to eavesdrop quantum-encrypted infor- 
mation. What makes the problem of gate programma- 
bility non trivial is that exact universal programmability 
of channels is impossible, as a consequence of a no-go 
theorem for programmability of unitary transformations 
by Nielsen and Chuang 3 . A similar situation occurs 
for universal programmability of POVM's 0J|S]- In both 
cases, it is still possible to achieve programmability prob- 
abilistically E3E1E], or even deterministically jjjjj, though 
within some accuracy. In establishing the theoretical lim- 
its to state-programmability of channels or POVM's the 
starting problem is to find the joint system-ancilla uni- 
tary or observable, respectively, which achieves the best 
accuracy for fixed dimension of the ancilla: this is ex- 
actly the problem that is addressed in the present paper. 
This problem turned out to be hard, even for low dimen- 
sion. Here we will give a solution for the optimal device 
for programming unitary channels for dimension two for 
both system and ancilla. On the other hand, as regards 
programming observables, we will give an upper bound 
for the optimal ancilla dimension d(e^ 1 ) versus the ac- 
curacy e—1 for programmable detectors. As we will see, 
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it turns out [S] that a dimension increasing poly- 

nomially with precision e~ Y is possible, and even a linear 
dependence is achievable for specific cases. This should 
be compared with the preliminary indications of an ex- 
ponential growth of Ref. However, even the linear 
dependence is still suboptimal. 

2 Statement of the problems 

Programmable unitaries We want to program uni- 
tary channels by a fixed device as follows 



\ t<T (p) = Tr 2 [V(p^a)V% 



(1) 



with the system in the state p interacting with an an- 
cilla in the state a via the unitary operator V of the 
programmable device (the state of the ancilla is the pro- 
gram). For fixed V the above map can be regarded as a 
linear map from the convex set of the ancilla states A to 
the convex set of channels for the system C. We will de- 
note by Vy.A the image of the ancilla states A under such 
linear map: these are the programmable channels. Ac- 
cording to the well known no-go theorem by Nielsen and 
Chuang it is impossible to program all unitary channels 
on the system with a single V and a finite-dimensional an- 
cilla, namely the image convex Vv,A C C is a proper sub- 
set of the whole convex hi of unitary channels and their 
convex combinations. This opens the following problem: 

Problem: For given dimension of the ancilla, find 
the unitary operators V that are the most efficient 
in programming unitary channels, namely which 
minimize the largest distance e(V) of each channel 
£ U from the programmable set T J v,A : 



e(V) = max min 8(tf, 3 s ) 



max min 6 iff, #y„). 

"U&AadA 

(2) 



As a definition of distance it would be most appropriate 
to use the CB-norm distance \€ — &\ CB - However, this 
leads to a very hard problem. We will use instead the 
following distance 



9) = y/1 - F{&, 0>), 



(3) 



where F( c (o , £P) denotes the Raginsky fidelity ^Uj , which 
for unitary map ^ = a i/ = U ■ U' is equivalent to the 
channel fidelity 

W^) = l^|Tr[G/C/]| 2 , (4) 

i 

where ^ = J^-Cj • C\. Such fidelity is also re- 
lated to the input-output fidelity averaged over all pure 
states F io (^,^>), by the formula F io (^, 9>) = [1 + 
dF(W, £?)]/{<! + 1). Therefore, our optimal unitary V 
will maximize the fidelity 

F(V) = min F(U, V), F(U, V) = maxF(W, &> Va ) 

(5) 

Programmable detectors The POVM of a measur- 
ing apparatus is a set of positive operators Pi ^ 0, i — 
1, . . . n, n < oo normalized to the identity X)i=i P = ^> 
which gives the probability distribution of the outcomes 
for each input state p via the Born rule 

p{i\p) ±Tt\pPi\. (6) 

Clearly, the most general programmable detector is de- 
scribed by an observable F = {Fi} jointly measured on 
system and ancilla. The probability distribution of the 
outcomes is given by 

p (T (i\p)=Tr[(p®a)F i \, Vi.Vp. (7) 

By taking the partial trace in Eq. J2J over the ancilla 
and using the polarization identity (Eq. J7J holds for all 
states) one obtains the POVM 

P„,i=Tc 2 [(I®<T)Fi]. (8) 

From Eq. J3J it follows that the convex set of states si 
of the ancilla is in correspondence via the map Pf.o- = 
Tr 2 [(J ® er)F] with a convex subset £ of 

the convex set of the system POVM's with n out- 
comes. The symbol ^py denotes the convex set of pro- 
grammable POVM's that can be achieved with fixed ob- 
servable F and varying state a £ si ' . The no-go theorem 
proved in E] states that for any fixed observable F 
the programmable set Pp.rf is strictly contained in &> n , 
since even just the observables cannot be programmed 
with a fixed observable F and a finite dimensional ancilla. 
We now restrict attention to programmability of observ- 
ables only, whence n = dim(Jff), the case of nonorthogo- 
nal POVM's simply resorting to program observables on 
a larger Hilbert space. In the following we will denote 
by O n the set of observables. The problem of measure- 
ment programmability can then be stated in mathemat- 
ical terms as follows 

Problem: For given dimension of the ancilla, find 
the joint observables F that are the most efficient 
in programming system observables, namely which 
minimize the largest distance e(F) of each observ- 
able P £ O n from the programmable set Pf.o--' 

e(F) = max min <5fP,Q) = max min <5(P. Pf a)- 
Peo„ QeP F a Peo„ a&A 

(9) 



We define the distance between two POVM's as the dis- 
tance between their respective probabilities, maximized 
over all possible states, namely 

<J(P,Q) =maxV|Tr[p(P l -Q l )]l- (10) 
p ~ 

The distance defined in Eq. (|10|) is hard to handle ana- 
lytically, whence we bound it as follows 

*(P,Q)<X)||Pi-Qi| <El P '-*l2> ( n ) 

i i 

where \A\ is the usual operator norm of A, and \\A\\ 2 = 
y^rpfL4] is the Frobenius norm. 

3 Programming qubit unitaries 

By some lengthy calculation we can obtain the Kraus 
operators for the map ^v,a(p) 

^V,cr(p) = J]] C nm pC nm , 

nm (12) 

C nm = J £e ie «* k \v* n )(v* m \*l^r n 

k 

where \v n ) denotes the eigenvector of a corresponding to 
the eigenvalue A„ and * denotes complex conjugation in 
the same fixed basis for which the operator have the 
same matrix elements as the matrix of coefficients of the 
eigenvector of V corresponding to eigenvalue e l9k . We 
then obtain 

I Tr[Cl m U}\ 2 =£V(«»-*0 Tr^t/tv^T^t^] 

nm kh 

= Tr[a T S(U,V)^S(U,V)} 

(13) 

where 

S(U,V) = Y,^ l6k ^lu^k- (14) 

fe 

and T denotes transposition in the canonical basis. The 
fidelity (0) can then be rewritten as follows 

F(U,V) = ±\\S(U,V)f. (15) 

The operator S(U, V) in Eq. (|14(l can be written as fol- 
lows 

S(U, V) = Tn[(J7 T ®I)V*\ . (16) 

Changing V by local unitary operators transforms 
S(U, V) in the following fashion 

S(U, (Wi <g> W 2 )V(W 3 ® Wi)) = W£S(W}UW£,V)WZ, 

namely the local unitaries do not change the minimum 
fidelity, since the unitaries on the ancilla just imply a dif- 
ferent program state, whereas the unitaries on the system 
just imply that the minimum fidelity is achieved for a dif- 
ferent unitary — say W^UW^ instead of U. 

For system and ancilla both two-dimensional, one can 
parameterize all possible joint unitary operators as fol- 
lows [TTlini (Wi <g> W 2 )V(W 3 <g> Wi), where 

V = cxp[i(aiC7i(g>cri T +a 2 o'2®o'2 1 + a 3 a 3 ®a 3 T )] . (18) 



The problem is now reduced to study only joint unitary 
operators of the form of Eq. I|18|) . It can be proved that 
the coefficients of its eigenvectors are the matrix elements 
of Pauli matrices <7j, j = 0, 1, 2, 3 where <7o = I, a\ = a x , 
o~2 = Cjn "3 = <y z . This means that we can rewrite 
S(U, V) in Eq. JTU) as follows 



with 



S{U,V)= l -Y.e 

3=0 



ai + a 2 + a 3 , 



2a, - 6n 



(19) 



(20) 



Through the derivation described in Appendix 0] we ob- 
tain that the fidelity minimized over all unitaries is given 
by 

1 



F(V) = -= min \t 
d ? 



where 



f 



3 
3=0 



(21) 



(22) 



tj = 



to, 1 < j < 3. 



The optimal unitary V is now obtained by maximiz- 
ing F(V). We need then to consider the decomposition 
Eq. I)18[l. and then to maximize the minimum among 
the four eigenvalues of S(U, V)^S(U, V). Notice that 
tj = V t HjiM e where H is the Hadamard matrix 



H = - 
2 



V 



-1 
-1 

1/ 



(23) 



which is unitary, and consequently X^I^jI 2 = 
J2j \e i93 \ 2 = 4. This implies that minj \tj\ < 1. We now 
provide a choice of phases 9j such that \tj\ = 1 for all j, 
achieving the maximum fidelity allowed. For instance, we 
can take 8 = 0, 9\ — tt/2, 9 2 — n, 63 — tt/2, correspond- 
ing to the eigenvalues — l,i for V. Another solution 
is 8q =0,6i = —ir /2,62 = 7T; ^3 = — tt/2. Also one can 
set 0, -> -6i. The eigenvalues of S(U, V)^S(U, V) are 
then 1,1,1,1, while for the fidelity we have 



F= max F(V) = i = 7, 
and the corresponding optimal V has the form 



V = cxp 



1 <J X ± <7, 



o"z) 



(24) 



(25) 



A possible circuit scheme for the optimal V is given in 
Fig. □ 

Such fidelity cannot be achieved by any V of the 
controlled-unitary form 



V 



2 

J2 V k®\^k)(*Pk\, 
fc=l 



<V>i|^2>=0, (26) 
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Figure 1: Quantum circuit scheme for the optimal uni- 
tary operator V in Eq. (|2*4"|) . VF Q denotes e*% aw . For 
the derivation of the circuit consider that <j x ® o~ x — 
C(cr x ® I)C and a z ® <j z = C(J ® cr x )C, where C de- 
notes the controlled- not. 



where Vi, V2 are unitaries on ~ C 2 . In fact, it is easy 
to see that in this case the fidelity is given by 

F(U, V) = ~| Tr[^£/]| 2 , ft = argmax| Tr^C/]!, 

(27) 

and for any couple of unitaries {14} there always exists 
a unitary U orthogonal to both {14}, whence F(V) = 
min Ue w(H) F(U,V) = 0. 

4 Upper bound on optimal size for pro- 
grammable detectors 

We will now derive an upper bound for the function 
d = d(e^ 1 ), where e = minp£(F), that gives the mini- 
mal needed dimension of the ancilla to achieve accuracy 
e _1 in programming observables for a finite-dimensional 
quantum system. Clearly the function d — d(e~ 1 ) must 
be increasing, since the higher is the accuracy e _1 , the 
larger the minimal dimension d needed for the ancilla, 
namely the "size" of the programmable detector. 

Consider now a <i-dimensional ancilla and a system- 
ancilla interaction U of the following controlled-unitary 
form 



u = J2w k ®\4> k ){<t>k\, 



(28) 



fe=i 



where {4>k} is an orthonormal complete set of vectors for 
the ancilla and W% are generic unitary operators on ffl . 
Consider then a POVM E = UFW of the form 



E t = \il}i){ipi\®lA, 



(29) 



where Ia denotes the identity operator on the ancilla 
space, and {V>fc} is a complete orthonormal set for the 
system. The observable to be approximated will then be 
written as follows 



P i = W*\il>i)(il> i \W, 



(30) 



W being a unitary operator on J^f, and we will scan 
all possible observables by varying W. For the program 
state of the ancilla we use one of the states (j>k, which 
give the POVM's 

Qi=Wl\T/j i )(ip i \W k . (31) 
This special form simplifies the calculation of the bound 



in Eq. Ijllfl . which becomes 



i(P,Q)<^^2(l-|{^* 



(32) 



and using the Jensen's inequality for the square root func- 
tion we have 



S(P,Q) < V2n\\W-W k \\ 



(33) 



Now we can always take d sufficiently large such that we 
can choose the d operators {VFfe} in the unitary trans- 
formation U in Eq. (|28(l in such a way that for each 
given W there is always a unitary operator Wk in the 
set for which \/2n \\W — Wk\ 2 is bounded by e. This will 
guarantee that for the given observable P correspond- 
ing to W there is a program state for the ancilla such 
that the POVM Q achieved by the programmable detec- 
tor is close to the desired P less than e. The set of all 
possible unitary operators W is a compact manifold of 



dimension h = 



We now consider a covering of 



the manifold with balls of radius 



centered at 



the operators Wk- This guarantees that any W would 
be within a distance — |^ from an operator Wk, which 
in turns implies that the accuracy of the programmable 
device is bounded by e via Eq. (|33|l . Using the volume 

h i_ 

of the ^.-dimensional sphere of radius r, we 



v ~ TTp+i) 

obtain the number of balls needed for the covering (for 
sufficiently small e, corresponding to the upper bound for 
the minimal dimension of the ancilla 



d ^ k(ji) 



i(n-l) 



(34) 



where n(n) is a constant that depends on n. Eq. (|34|l 
gives an upper bound for the dimension d which is poly- 
nomial versus the accuracy 

For qubits, the observable has only two elements, Pq = 
\fp)(tp\ and Pi = \tp±)(ip± \ = I — Pq, and the distance in 
Eq. H1(J[1 can be evaluated analytically as follows 



<J(P,Q) 



max2|Tr[p(Po 

p 



Qo) 



(35) 



As regards now the programmability of all POVM's 
(i. e. including the nonorthogonal ones), just notice 
that one just needs to be able to program only the ex- 
tremal POVM's in since their convex combinations 
will corresponds to mixing the program state or to ran- 
domly choosing among different detectors. Then, since 
their maximum number of outcomes is n 2 , the extremal 
POVM's have Naimark's extension to observables in di- 
mension n 2 , whence we are reduced to the case of pro- 
grammability of observables in dimension n . 

We will now give a programmable detector for qubits 
that achieves an accuracy that is linear in d. For the an- 
cilla we use a generic d-dimensional quantum system, and 
relabel the dimension in the angular momentum fashion 
d = 2j + 1. The idea is now to design a programmable 
detector in which the unitary transformation correspond- 
ing to the observable {Pi} in Eq. i(3"U|) is programmed 



by covariantly changing the program state of the ancilla. 
By labeling unitary transformations by a group element 
g € §U(2), we write the observable to be programmed 
as P = K?||)(f \Vg wnere {Vg} = (h) 18 a unitary irre- 
ducible representation of SU(2) with angular momentum 
|, whereas the program state will be written as W g crW£, 
with {VF 9 } = (j) a unitary irreducible representation of 
§U(2) on the ancilla space with angular momentum j. 
As already noticed, without loss of generality we can al- 
ways choose the state a as pure. We will now show that 
a good choice for the program state is a = 
{ \j, m) } denoting an orthonormal basis of eigenstates of 
J z in the irreducible representation with angular momen- 
tum j. The tensor representation {V g <g> W g } = | j can 
be decomposed into the direct sum of two irreducible 
representations \ ® j = j+ j_, where j± = j ± F° r 
the POVM F of the programmable detector we will use 
Fq = Z + and Ft = Z- , Z± denoting the orthogonal pro- 
jector on the invariant space for angular momentum j± 



Fn = 



3+ 



(36) 



-3+ 



Using the invariance (V g <g> W g )F (V<j ® W%) = F , we 
can write the programmed POVM as follows 

Qo =Tr A [(I® W%\jJ)(j,j\W g )F ] 

=V2Tr A [(I®\j,j){j,j\)F ]V g 

-v (\i i I _i i I i _ 

v 3 y\ 2 ' 2 / \ 2 ' 2 I ~ 2j+l I 2 ' 2/\2' 



(37) 



' a 7 



where we used the only non vanishing Clebsch- 
Gordan coefficients j}\ 2 = lj and 

|0+,i_|i,-i)|i,i)| 2 = 5^3-. Clearly, Q - P Q = 
2j+i^ r g\l' — D(l' — il^s' whence according to Eq. l|3"K|l 
the accuracy is given by S(P, Q) = 2/d. The scaling of 
the dimension with the accuracy is then linear 

d = 2E~\ (38) 

whereas the bound (|3~4"jl would be quadratic d oc e -2 . 
Sublinear growth of d versus e" 1 is not excluded in gen- 
eral, but is not possible for the present model. 

Appendix: Derivation of the minimum fi- 
delity for the unitary V 

Starting from Eq. JTHl we will obtain Eq. pi)l. The 

unitary U belongs to SU(2), and can be written in the 
Bloch form 

U = uqI + in ■ a , (39) 
with rife € K and n\ + |n| 2 = 1. Using the identity 

(JjUlCTj = Cjiai, 6j0 = 6jj = 1, 6jl = — 1,1 ^ 0,j, 

(40) 

we can rewrite 



where 



S(U, V) = h I + ri- ct, 



1 3 

--tjUj, 0<j<3, t = -'^2> 



tj =e 



3=0 

to, 1 < j < 3, < j < 3, 



(41) 



(42) 



and we will use the exponential representation for the 
complex number tj = \tj\e lc ^ j . It is now easy to evaluate 
the operator S(U, V)^S(U, V). One has 



S(U,V)^S(U,V) =v Q I + v 



a, 



\no\ 2 



23(n n*) + n* x n 



(43) 



Now, the maximum eigenvalue of S(U, Vy S(U, V) is i>o + 
\v\, and one has 



\v\ 2 = 2 In/I^fsin 2 ^- 

i,j=0 

whence the norm of S(U, V) is given by 



\\S(U,V)f = u-t- 



Tu . 



(44) 



(45) 



where u 
and Tij = \U 



tQ, 11,1, , 
\ti 



' sm 



2 * "3 

2 / j 



(|*o| a , l*i| a , Ital 2 , ]t 3 | 2 ), 
Notice that the uni- 



tary U which is programmed with minimum fidelity in 
general will not not be unique, since the expression for 
the fidelity depends on {nf}. Notice also that using the 
decomposition in Eq. (|18|l the minimum fidelity just de- 
pends on the phases {Oj}, and the local unitaries will 
appear only in the definitions of the optimal program 
state and of the worstly approximated unitary. One has 
the following bound on the expression in Eq. 14511 



t + V u ■ Tu > u ■ t > min \tj 



(46) 
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and the bound is achieved on one of the four extremal 
points ui = 5ij of the domain of u which is the convex 
set {u, uj > 0, J2j u j = 1} (the positive octant of the 
unit four dimensional ball Si.). This proves the content 
of Eq. (JUJ. 
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